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ABSTRACT 

A  model  is  developed  which  relates  the  structure  of  a 
combat  unit  to  its  expected  effectiveness  during  a  time  inter- 
val.  Sensitivity  to  variations  of  structure  is  incorporated 
into  the  model  by  adapting  effectiveness  concepts  from  reli- 
ability theory,  and  by  constructing  the  attrition  process  as 
a  continuous  parameter  Markov  chain.   The  model  is  highly 
sensitive  to  the  judgment  of  the  user,  due  to  parametric 
flexibility  both  in  the  attrition  process,  and  effectiveness 
measurement.   Two  solution  methods  are  given.   Examples  are 
employed  to  illustrate  the  model's  use,  and  proposals  for 
refinement  are  offered. 
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I.   INTRODUCTION 

The  purpose  of  any  operational  model  is  to  aid  a  decision- 
maker by  relating  his  possible  alternatives  to  their  consequent 
results,  thus  increasing  the  knowledge  upon  which  his  decision 
is  based.   Therefore  a  model,  in  addition  to  validly  represent- 
ing the  process  or  system  involved,  must  be  capable  of  dis- 
playing the  contrast,  in  terms  of  meaningful  criteria,  of  even 
subtle  variations  of  choice.   Sensitivity  to  variables  outside 
the  area  of  the  current  decision  is  not  particularly  valuable, 
and  in  fact  may  produce  detrimental  complication.   In  short,  a 
model  should  ideally  be  a  reflection  of  the  problem  of  its  user 
It  is  for  this  reason  that  distinct  decision  problems  often 
require  the  use  of  different  models  of  a  single  system. 

This  dependence  of  model  upon  problem  is  nowhere  more 
apparent  than  in  the  study  of  combat  dynamics.   As  an  example 
a  single  battle  situation  may  be  modeled  using  game  theory  if 
the  decision  involves  strategy,  or  using  search  theory  if 
alternative  fire  distributions  are  being  studied,  or  using 
Lanchester  theory  if  weapon  characteristics  or  tactical  move- 
ments are  being  considered.   Each  modeling  approach  has  been 
found  useful  for  certain  problem  types  but  useless  for  others. 

The  purpose  of  this  thesis  is  to  propose  an  alternate 
modeling  viewpoint  for  those  problems  in  which  the  structure 
of  a  combat  unit  must  be  related  to  its  effectiveness.   The 
term  structure  is  meant  to  denote  the  personnel  and  equipment 


of  the  unit,  the  roles  of  these  elements,  as  well  as  the 
channels  employed  to  receive  and  transmit  information  and 
instructions.   As  an  example,  a  decision-maker  who  must 
choose  for  field  test  a  small  number  of  prototype  tank  killer 
teams  from  all  the  infinite  organizational  possibilities, 
faces  a  structural  decision  problem,  and  might  profit  from 
the  model  to  be  proposed. 

The  importance  of  structure  derives  from  the  fact  that 
when  a  leader  or  communicator  becomes  a  casualty,  his  removal 
from  the  structure  produces  twin  transformations  in  both  the 
effectiveness  and  vulnerability  of  the  unit.   These  transforma- 
tions are  determined  at  least  in  part  by  the  casualty's  posi- 
tion in  the  organization.   As  an  illustration  consider  a  rifle 
squad  in  defense  whose  squad  leader  has  just  become  the  first 
casualty.   His  function  is  assumed  by  a  fire  team  leader,  who 
must  now  receive  and  transmit  orders,  direct  the  fire  of  his 
squad  and  search  his  sector  for  targets.   The  effectiveness  of 
the  squad  is  reduced  not  only  because  one  fewer  weapon  is 
firing  but  because  a  less  highly  trained  leader  is  in  charge. 
The  vulnerability  has  increased  for  the  team  leader  because  he 
now  must  increase  his  exposure  in  order  to  direct  the  whole 
squad  and  observe  its  whole  sector.   The  vulnerability  of  the 
whole  squad  is  perhaps  increased  because  it  is  being  directed 
with  less  expertise. 

Thus,  to  know  the  condition  of  a  unit  it  is  not  enough  to 
know  the  number  of  casualties  without  knowing  who  they  are. 
Perhaps  the  strongest  testimony  to  this  viewpoint  is  the 


universal  military  principle  that  enemy  leaders  and  radio 
operators  are  the  targets  of  highest  priority  in  any  type  of 
engagement,  and  the  simultaneous  effort  to  conceal  visual 
signs  of  their  identity.   If  the  validity  of  this  approach 
has  been  justified  the  next  step  is  to  outline  the  model 
requirements . 

The  decision-maker  will  be  assumed  to  have  available  to 
him  a  facility  for  military  judgment  which  will  include  his 
own  assessment  of  unit  and  individual  vulnerability,  and  his 
own  criteria  for  measuring  unit  effectiveness  on  the  basis  of 
its  surviving  structure.   The  mission  of  the  model  will  be  to 

1)  provide  a  framework  in  which  the  decision-maker's 
assessment  of  threat  can  be  specifically  quantified, 

2)  translate  the  threat  assessment  into  a  stochastic 
description  of  the  behavior  of  the  candidate 
organization, 

3)  provide  a  framework  for  expressing  the  decision- 
maker's effectiveness  criteria  as  an  algebraic 
function  of  the  unit  condition/ 

4)  derive  a  stochastic  description  of  effectiveness  of 
the  candidate  organization,  using  the  effectiveness 
function. 

Fortunately,  the  complete  process  of  model  construction  is 
unnecessary,  since  there  already  exists  an  established  field 
of  study  which  treats  the  concept  of  structure.   The  field  is 
stochastic  reliability  theory,  which  has  historically  devel- 
oped in  answer  to  the  need  for  relating  the  design  of  complex 
systems  of  mechanical  and  electronic  components  to  system 
effectiveness.   The  comparison  of  a  system  of  components  to 


a  unit  of  fighting  men  may  seem  rash.   However,  through  a 
process  of  adaptation,  a  model  will  be  derived  which  may  be 
considered  more  consistent  while  still  taking  advantage  of 
reliability  concepts. 

Additional  reliance  on  the  theory  of  continuous  parameter 
Markov  processes  will  permit  treatment  of  significant  inter- 
dependencies  in  the  attrition  process. 

The  resultant  model  is  one  with  capacity  for  immediate 
application,  as  well  as  promising  opportunities  for  improve- 
ment and  refinement. 


II.   RELIABILITY  CONCEPTS 

Reliability  theory  deals  with  collections  of  inter- 
dependent components.   These  collections  are  called  systems. 
A  component  is  any  entity  which  can  be  described  at  any  time 
as  either  functioning  or  failed,  depending  upon  its  ability 
to  perform  its  assigned  mission.   Components  are  generally 
assumed  to  be  in  one  of  these  two  conditions.   No  partially 
failed  conditions  are  considered. 

The  system  which  is  discussed  in  this  thesis  is  an  engaged 
combat  unit  whose  components  are  included  in  one  or  more  of 
the  following  categories: 

1)  personnel, 

2)  weapons  and  equipment, 

3)  channels  of  information  and  communication. 

The  number  of  components  in  the  unit  is  represented  by 
the  variable  n  called  the  order  of  the  unit.   Each  component 
is  identified  with  an  integer  in  the  set  {l,...,n},  called  its 
index . 

The  performance  random  variable  of  the  ith  component  is 

defined  as  follows: 

X.  (t)  =  0    if  the  ith  component  is  in  the 
failed  condition  at  time  t, 

=1    if  the  ith  component  is  in  the 

functioning  condition  at  time  t. 

The  component  reliability  of  the  ith  component  is  defined 

as  follows: 


ri(t)  =  P[Xi(t)  =  1], 

Since  there  are  n  components  of  interest,  it  is  an  obvious 
extension  to  define  the  performance  random  vector  as  follows: 

X(t)  =  (x1(t) ,.. .,xn(t))  . 

The  outcome  space  of  X(t)  is  the  set  of  all  binary  n-vectors. 

Call  this  set  ft  .   This  set  has  N  members  where  N  =  2  . 

n 

It  is  often  necessary  to  partition  the  outcome  space  ft 

into  m  subsets,  a, .....a  ,  where  1  <  m  <  N,  and  define  a 

1      m  —   — 

partition  function  by  assigning  a  unique  real  number  b.  to 
all  vectors  in  each  subset  a. ,  so  that 

a  (x)  =  b .     if  and  only  if    x  e  a. . 
—     i  J       —    i 

Each  vector  is  a  member  of  exactly  one  of  the  subsets. 

The  first  partition  function  of  interest  is  the  state 

function  s(x),  which  is  constructed  by  partitioning  ft   into 

N  subsets,  each  containing  one  vector.   The  value  assigned  to 

each  one-member  subset  is  one  of  the  integers  1,...,N. 

Therefore 

s (x)  =  j     if  and  only  if     x  e  a.  and  b.  =  j . 

It  is  assumed  here  that  in  all  cases 

s (1)  =  1     and 

s  (0)  =  N    where 

1  =  (1, . . . , 1)     and 

0  =  (0, ... ,0)  . 
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The  composition  of  the  state  function  with  the  perfor- 
mance random  vector  results  in  definition  of  the  state  random 
variable 

S(t)  =  s(X(t))  . 

The  common  practice  in  reliability  theory  is  to  consider 
the  entire  system  to  be  either  in  the  functioning  or  failed 
condition  at  any  time.   Predictably  the  system's  condition  is 
described  by  a  binary  random  variable.   However,  when  the  sys- 
tem is  a  combat  unit,  such  a  measure  is  too  coarse  to  be  of 
value.   Still,  since  this  binary  concept  is  of  value  in  the 
adaptation  process,  it  is  introduced  below. 

The  structure  function  0(x)    of  a  system  of  order  n  is  a 
partition  function  constructed  by  forming  two  subsets: 

1)  a, ,  the  set  of  all  vectors  x  in  ft  ,  such  that  if 

1'  —     n' 

X(t)  =  x,  then  the  system  is  functioning, 

2)  a~,  the  set  of  all  vectors  y_  in  9,    ,  such  that  if 
X(t)  =  y,  the  system  is  in  the  failed  condition. 

Consistent  with  the  use  of  binary  indicators, 

p'(x)  =  l(=b,)     if  and  only  if     x  z    a, 
=  0(=b„)     if  and  only  if     x  z    a~. 
Composing  the  structure  function  with  the  performance 
random  vector  results  in  definition  of  the  structure  random 
variable 

*(t)  =  0(X(t))  . 
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The  system  reliability  is  defined  to  be 

R(t)  =  P[$(t)  =  1] . 

If  more  than  one  system  is  being  considered  then  it  is  neces- 
sary to  place  a  subscript  on  the  structure  function,  struc- 
ture random  variable  and  system  reliability. 

A  system  is  called  coherent  if  its  structure  function 
satisfies  three  conditions: 

1)  jz*(_0)  =  0, 

2)  0(1)  =  1, 

3)  If  x  <  y_,   then   p'(x)  £  jzJ  (y_)  . 

Essentially  a  coherent  system  is  one  which  fails  if  all  its 
components  fail,  functions  if  all  its  components  function,  and 
whose  condition  is  not  impaired  by  improving  the  condition  of 
any  of  its  components. 

A  component  is  called  irrelevant  if  its  condition  never 
has  any  effect  on  the  value  of  the  structure  function. 

Components  are  assumed  to  be  of  two  varieties.   Renewable 
components  alternate  between  the  functioning  and  failed  condi- 
tions.  For  example,  tactical  radios  are  communication  channels 
whose  functioning  is  interrupted  by  atmospheric  disturbance  or 
jamming.   Nonrenewable  components  never  return  to  the  function- 
ing condition  after  entering  the  failed  condition.   Personnel 
are  nonrenewable  components  of  a  unit. 

A  common  graphic  method  of  illustrating  the  relationship 
of  component  performance  to  system  performance  is  the  use  of 
the  block  diagram.   Examples  of  these  diagrams  are  shown  in 
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Figures  1-4.   The  method  is  derived  from  the  common  practice 
of  representing  electrical  circuits.   Each  block  corresponds 
to  the  like  indexed  component  of  the  system.   A  potential 
flow  is  imagined,  attempting  to  pass  through  the  diagram 
along  the  connecting  arcs.   The  flow  is  able  to  pass  through 
a  component  block  only  if  the  component  is  functioning.   If 
the  flow  is  able  to  traverse  the  diagram  over  any  path  then 
the  system  is  functioning  and  its  structure  function  equals 
one.   Irrelevant  components  are  represented  by  broken  blocks 
totally  unconnected  to  the  system. 

Figure  1  shows  the  block  diagram  for  the  system  of  three 
components  which  functions  if  and  only  if  all  three  components 
are  functioning.   This  is  called  a  series  system.   Figure  2 
shows  a  system  which  fails  if  and  only  if  all  its  components 
fail.   This  is  called  a  parallel  system.   Figure  3  shows  a  sys- 
tem which  functions  if  and  only  if  at  least  two  of  the  compo- 
nents are  functioning.   This  is  called  a  two-out-of-three 
system.   As  in  this  case,  a  single  component  may  be  represented 
by  more  than  one  block,  with  the  understanding  that  all  blocks 
with  the  same  index  are  always  in  the  same  condition.   Figure 
4  shows  a  system  with  components  one  and  two  in  parallel,  with 
component  three  irrelevant.   It  should  be  remembered  that  the 
block  diagram  shows  only  structural  relationship,  and  implies 
nothing  concerning  actual  location. 

A  final  concept,  not  unique  to  reliability,  but  useful 
here,  is  that  of  Laplace  transformation.   Given  a  function 
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Figure  1.   Series  System 


14 


1 

2 

3 

i 

Figure  2.   Parallel  System 
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Figure  3.   Two-Out-of-Three  System 
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Figure  4.   Parallel  System  with  Irrelevant  Component 
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c(t),  the  Laplace  transform  of  c(t)  is  defined  as  follows: 

oo 

L  (c(t)}  =  c(s)  =   /  exp(-st)  c(t)  dt, 

o 

a  function  of  s.   The  function  c(t)  is  called  the  inverse  of 
the  Laplace  transform  c(s). 
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III.   MEASUREMENT  OF  EFFECTIVENESS 

In  order  to  make  any  description  of  unit  performance 
useful  to  a  decision-maker,  it  is  necessary  to  translate  that 
description  into  some  measure  of  effectiveness  or  goal  attain- 
ment.  Such  a  measure  is  useful  only  if  it  is  consistent  with 
the  decision-maker's  standards  of  value,  and  is  in  a  form 
sufficiently  simple  to  permit  reasonable  comparison  of  compet- 
ing alternatives.   The  method  developed  here,  of  constructing 
such  a  measure  begins  with  the  designation  by  the  user  of  a 
relatively  simple  instantaneous  indicator  function  based  upon 
his  intuitive  value  judgments.   This  function  then  serves  as 
a  basis  for  derivation  of  a  one  dimensional  deterministic 
function,  which  measures  unit  quality  over  a  variable  time 
interval . 

The  actual  effectiveness  of  an  engaged  unit  at  a  point  in 
time  is  a  very  complex  function  of  numerous  variables,  includ- 
ing friendly  and  enemy  force  levels,  and  many  immeasurable 
physical  and  psychological  factors.   The  simplifying  assumption 
is  made  here  that  effectiveness  depends  only  upon  the  vector  of 
component  performance  and  the  length  of  time  elapsed  since  the 
beginning  of  the  engagement.   The  enemy  performance  is  not 
considered  directly,  but  the  designation  of  an  effectiveness 
function  should  reflect  consideration  of  a  hostile  enemy  envi- 
ronment, based  on  a  previously  determined  scenario.   The 
friendly  unit  is  always  assumed  to  have  all  components  func- 
tioning at  time  zero. 
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The  decision-maker's  first  step  is  to  define  a  partition 

function  on  ft  .   The  set  is  partitioned  into  subsets  a., 
n  1 ' 

i=l,...,m/  where  vectors  in  the  same  subset  are  judged  to 
represent  equal  effectiveness  unit  conditions.   Vectors  in 
different  subsets  represent  outcomes  which  are  not  equally 
effective.   This  partition  is  performed  without  reference  to 
time.   To  each  subset  a.  is  assigned  a  dimensionless  real 
number  g.,  called  its  absolute  effectiveness  level.   The 
function  thus  defined  is  called  the  absolute  effectiveness 
function  g(x)  where 

g(x)  =  g.     if  and  only  if     x  e  a. , 

and  the  partition  is  constructed  so  that 

*1    <  g2  K    •••  <  gm ' 
It  is  assumed  without  serious  limitation  that 

(3-1)  g(0)  =  0, 

(3-2)  g(l)  -  gm, 

(3-3)       if   x^y     then     g(x)  <  g  (y_)  . 

These  assumptions  merely  require  that  a  completely  destroyed 

unit  be  considered  completely  ineffective,  that  a  completely 

undamaged  unit  be  considered  maximally  effective,  and  that 

loss  of  a  component  never  improve  effectiveness. 

Since  the  vector  members  of  ft   and  the  integers  in  the 

n 

set  {1,...,N}  are  placed  in  one-to-one  correspondence  by  the 
state  function  s  (x)  ,  it  is  possible  to  designate  a  companion 
function  to  g(x)  called  the  alternate  effectiveness  function 
defined  on  {1,...,N}  such  that 
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d(j)  =  g.     if  and  only  if     s (x)  =  j   and   x  e  a.. 
It  follows  that 

g(x)  =  d(s(x))  . 

This  function  is  useful  in  later  calculations. 

The  construction  of  the  actual  effectiveness  function  is 
of  extreme  importance  to  the  successful  employment  of  the 
model.   The  user  should  take  special  care  to  assure  that  the 
absolute  effectiveness  levels  truly  represent  his  judgment, 
not  only  with  respect  to  ordinal  relationships,  but  to  ratio 
relationships  as  well. 

A  unit's  real  effectiveness  depends  upon  how  long  it  has 
been  engaged.   One  factor  to  be  considered  is  the  change  in 
enemy  situation  due  to  losses  or  reinforcement  as  provided  in 
the  scenario.   Another  consideration  is  that  in  missions  of 
fixed  length,  such  as  delaying  actions,  casualties  are  less 
damaging,  the  later  they  occur.   To  allow  for  such  trends,  the 
user's  next  step  is  to  define  a  time  variance  function,  f (t) , 
and  a  resultant  actual  effectiveness  function, 

h(x,t)  =  g(x)  f  (t)  . 

If  in  a  specific  application  the  variance  effect  is  not  con- 
sidered, then 

f  (t)  =1,  t  >  0 

and  therefore 

h(x,  t)  =  g(x)  ,  t  >.  0. 
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The  function  of  g  (x)  can  be  expressed  in  another  equiva- 
lent manner  which  shows  the  close  relationship  with  reli- 
ability theory.   Esary,  Proschan,  and  Walkup  [Ref.  3]  point 
out  that  a  function  such  as  g(x),  having  a  finite  range  can 
be  written  as 

(3-4)    g(x)  =  (g2-g1)  ^(x)  +  (g3-g2)  02{-)   +  ••• 

+  (g  -g  i )  0     i  (x) 

^m  ^m-1   >m-l   — 
where  the  functions  0 . (x)  are  defined  such  that 

(3-5)     o\  (x)  =0     if  and  only  if     g(x)  <_   g. 

(3-6)     0. (x)  =1     if  and  only  if     g (x)  >  g. . 

It  is  now  shown  that  each  of  the  functions  6 .  (x)  is  the  struc- 

l  — 

ture  function  of  a  unique  binary  coherent  system  composed  of 
the  components  of  the  unit,  in  which  at  least  one  of  the 
components  is  not  irrelevant.'   Furthermore  the  systems  are 
related  in  such  a  way  that  no  system  is  functioning  if  any 
lower  indexed  system  has  failed,  and  no  system  is  failed  if 
any  higher  indexed  system  is  functioning. 
Theorem  1_ 

The  functions  0  . (x)  in  (3-4)  are  structure  functions  of 
coherent  systems  of  the  n  components  of  the  parent  system. 
Proof 

It  must  be  shown  that  for  i=l,...,m, 

1)  0i(O)  =  0, 

2)  0i(l)  =  1, 


3)   if  x  <_  y_    then    0\  (x)  <  o^  (y_)  . 
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The  proof  is  by  contradiction.   The  above  three  points  are 
each  in  turn  assumed  false,  and  a  result,  contradicting  one 
of  the  assumptions  made  in  defining  g (x) ,  is  shown  to  follow. 

1)  Assume  that  0.  (0)    ^  0    for  some  i.   Then  it  must  be 

that  0.(0)  =  1#  but  then  from  (3-4) 

g(o)  >  (g1+1-g±)  >  o 

contradicting  (3-1).   Thus  0.(0)    =  0,  for  i  =  l,...,m-l, 

2)  Assume  that  0. (1)    ^    1    for  some  i.   Then  it  must  be 
that  0.(1)  =  0,  but  then  from  (3-4) 

g(i)   <   (g2-gi)  +  ...  +  (gm-gm_!)  -  gm,   or 

*<!>  <  *m 
contradicting  (3-2).   Thus  0.(1)    =  1,  for  i=l,...,m-l. 

3)  Assume  that  there  exist  vectors  x  and  y  in  Q    ,  such 
that  x  <_  y_  and  that  for  some  i,  0.   (x)  >  0.   (y_)  .   Since 
the  function  is  binary,  0.  (x)  =  1,  and  0.   (y_)  =0.   By 
(3-6),  0 . (x)    =  1    implies  that    g(x)  >  g..   By  3-5), 

0.   (y_)  =  0    implies  that    g(y_)  <.  g  •  •   Therefore 

g  (y )  <  g(x_)     and    x  <_  y 

contradicting  (3-3)  .   Thus  if  x  ±   y ,  then 

0i(x)  <_  0i(y)  ,    for  i  =  1,...,  m-1.  || 

Theorem  2 

The  functions  0 . (x) ,    i  =  1,...,  m-1,  are  related  as  follows 
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1)  if  0,{x)    =   1,  then  2$  .  (x)  =  1,  for  j  =  1,...,  k-1, 

K  —  J     — 

2)  if  0,   (x)  =  0,  then  0 . (x)  =  0,  for  j  =  k  +  1,...,  m-1. 

K  —  J  — 

Proof 

The  proof  is  again  by  contradiction. 

1)  Let  PVCx)  =  1  for  some  vector  x.   If  k=l  then  the 

result  is  trivially  true.   If  k  >  1,  assume  that 

0 . (x)  =  0  for  some  j  <  k.   By  (3-6),  0,  (x)  -  1  implies 
3  —  K  — 

that  g(x)  >  g,  .   By  (3-5),  0 .  (x)  =  0  implies  that 

—       K  J     — 

g(x)  <.  g  .   Since  g.  <  g,  ,  it  follows  that 

—  3  j     K 

g(x)  <  gk  <  g(x)  , 
an  obvious  contradiction.   Therefore  0 . (x)  =  1  for 
j  =  1, ... ,k-l. 

2)  Let  0,  (x)  =  0  for  some  vector  x.   If  k  =  m-1,  the 

result  holds  trivially.   If  k  <  (m-1),  assume  that 

^ .  (x)  =  1  for  some  j  >  k.   By  (3-5),  o\  (x)  =  0  implies 
3  —  K  — 

that  g(x)  <_  gv .   By  (3-6)  ,  0.  (x)  =  1  implies  that 

g(x)  >  g..   Since  g.  >  g,  ,  it  follov:s  that 

—  3  3  K 

g(x)  <  g  .  <  g(x)  , 

an  obvious  contradiction.   Therefore  p'-(x)  =  0,  for 

j  =  k  +  1,  .  .  .  ,  m-1.  || 
Theorem  3_ 

For  any  vector  x,  g(x)  =  g.,  if  and  only  if  0 . (x)  =  0,  and 
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Proof 

1)  Let  g(x)  =  g^  for  some  x.   Then  by  (3-5),  g(x)  = 

gi  -  gi  imPlies  that  &^W    =  °-   BY  (3-6),  g(x)  = 
gi  >  gi-l  imPlies  that  ^i_i^x)  =  1. 

2)  Let  0i(x)  =  0,  and  ^i_-,(x)  =  1  for  some  x.   By  (3-5) 
0 .  (x)  =  0  implies  that  g(x)  <.  g.  .   By  (3-6), 
0i_1(x)  =  1  implies  that  g(x)  >  g.  ,.   Therefore  it 
follows  that 

gi-l  <  g(-)  -  gi' 
but  the  only  value  in  the  range  of  the  function  v/hich 
satisfies  the  inequality  is  g. .   Therefore  g(x)  =  g-.| 
The  ultimate  importance  of  this  method  of  expressing  g(x) 
is  that  the  study  of  a  very  complicated  multivalued  absolute 
effectiveness  function  can  be  simplified  to  the  study  of  a 
linear  combination  of  simple  dependent  binary  structure  func- 
tions of  the  type  commonly  treated  in  classic  reliability. 

It  follows  that  the  actual  effectiveness  function  can  be 
written  as  follows: 

h(x,t)  =  f(t)  {(g2-gi)  01(x)    +    ...  +  (g^g^i)  ^m_i^>}- 

Thus  far  the  functions  introduced  are  deterministic  with 

domain  fi    .   Since  Q      is  also  the  outcome  space  of  the  perfor- 
n  n  c  r- 

mance  random  vector,  X(t) ,  each  of  the  functions  can  be 
composed  with  this  vector  to  define  analogous  random  variables 
The  absolute  effectiveness  random  variable  is  defined  by  the 
relation 

G(t)  =  g(X(t)  )  , 
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or    equivalently 

G(t)    =    (g2-gi)    Vt)    +    ...    +    (g^g^)    Vl(t>< 
where 

$i(t)    =   0i(X{t))t    i   =   1,...,    m-1. 

The  actual  effectiveness  random  variable  is  defined  by  the 
relation 

H(t)  =  f  (t)G(t)  . 

These  effectiveness  random  variables  have  the  advantage 
that  they  are  defined  through  a  process  which  assures  their 
consistency  with  the  decision-maker's  standards.   Their  pri- 
mary limitation  is  that  they  represent  only  instantaneous 
effectiveness.   That  is  for  every  nonnegative  value  of  t 
there  is  defined  a  distinct  random  variable  H(t) .   In  Figure 
5,  the  bold  horizontal  line  depicts  a  possible  curve  of  out- 
comes of  H(t),  where  f(t)  is  assumed  to  equal  one  for  all  non- 
negative  t. 

Since  combat  engagements  are  extended  over  indefinite  time 

periods,  it  is  desirable  to  aggregate  the  effectiveness  curve 

into  a  variable  which  measures  accumulated  effectiveness  over 

an  interval.   To  fill  this  need  the  value  random  variable  is 

defined  by  the  relation 

t  t 

V(t)  =    /   H(s)  ds  =   J   G(s)  f(s)  ds. 
o  o 

In  Figure  5,  the  value  random  variable  is  represented  by  the 

shaded  area  below  the  curve  of  H(t). 
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Figure  5.   Example  Outcome  of  Actual  Effectiveness  Random 
Variable  and  Value  Random  Variable 
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The  value  random  variable  is  a  measure  of  the  type 
desired,  but  in  general  it  will  have  a  complicated  distribu- 
tion, too  unwieldy  to  be  of  practical  use  in  decision  making. 
A  relatively  simple  measure  of  effectiveness,  which  still 
preserves  many  advantages  of  the  value  random  variable,  is 
the  quality  function  defined  by  the  relation 

q(t)  =  E[V(t)], 
the  mean  or  expected  value  of  the  value  random  variable. 

To  compute  the  quality  function  in  a  specific  case,  note 

that 

t 
q(t)  =  E    I  H(u)  du 
L  o 


E[H(u)]  du 


/  f  (u)  E[G(u)  ]  du, 
o 

where  the  interchange  of  expectation  and  integration  in  the 
second  line  is  justified  by  application  of  Fubini ' s  Theorem, 
since  H (u)  is  a  family  of  nonnegative  random  variables.   The 
validity  of  this  interchange  is  established  by  Royden  [Ref. 
12]. 

There  are  two  approaches  for  calculating  q(t).   The  first 
arises  from  the  fact  that 
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[G(uJ]    =     £      g(x)    P[X(u)     =    x    |X(0)     =    1] 


X 

N 
=       £d(i)    p1±(u)f 
i=l 

where    . 

p1JL    (u)    =    P[S(u)    =    i    |     S(0)    =    1] 

is  called  the  interval  probability  function  from  state  one  to 

state  i.   Thus  the  first  result  is  that 


q(t)  =  I   f(u) 
o 


N 


)  d(i)  p,  . (u) 
L  i=l 


du 


d(i) 


N 

L 

i=l 
N 


f(u)  P1:i_(u)  du 


(3-7 


where 


)  =   ^d(i)  I(i't} 


i=l 


(3-8)   I(i,t)  =   /   f(u)  P1;L(u)  du. 
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The  second  approach  follows  from  relation  (3-4) , 

m-1 
E[G(u)j  =  E[  ^      ^j  +  i-^j)  *.(u)J 
j-l  J    . 

m-1 

=  E  (gj+l-gj}  BI»j^)3 
j-l 

m-1 

=  E   (gj+l-gj}  P[0J(U)  =  1] 
j-l 

m-1 

■  I  (*j+.rv  vu)' 

j-l 

where  R.(u)  is  the  system  reliability  of  the  jth  associated 
system,  computed  on  the  basis  of  all  components  functioning 
at  time  zero.   And  the  second  result  is  that 


r  m-1 


q(t)  =  J   f(u) 
o 

m-1 


E  (^i-v  vu) 


L  j=l 


du 


j-l 

m-1 


f (u)  R. (u)  du 


L-   O 


(3-9) 


where 


/  , 
j=l 


(gj+l"gj}  J(3/t) 


(3-10) 


J(j,t)  =   /  f(u)  R. (u)  du. 


Now  since  the  functions  f(t)  and  g(x)  or  d(i)  are  assumed 
already  defined  by  the  decision-maker,  all  that  is  required  for 
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calculation  of  the  quality  function  is  either  set  of 
functions  {p11(u),...f  Pj_n(u)K  or  (R^di)  , . .  . ,  Rm_1(u)^«   The 
determination  of  the  functions  is  treated  in  Section  IV. 

The  quality  function  is  the  primary  measure  of  effec- 
tiveness proposed  here.   It  is  a  one  dimensional  deterministic 
function,  yet  it  is  based  upon  consideration  of  the  random- 
ness of  unit  performance.   Since  it  is  derived  from  functions 
defined  by  the  decision-maker,  it  is  consistent  with  his 
implicit  judgment.   And  since  it  is  a  function  of  time,  it 
displays  sensitivity  to  the  length  of  engagement. 

In  order  to  illustrate  the  use  of  the  quality  function  in 
predicting  unit  effectiveness,  this  section  will  be  closed  with 
an  example. 
Example  1 

Consider  a  tank  killer  team  consisting  of  the  following 
components : 

COMPONENT  INDEX 

Team  Leader         1 
Gunner  2 

Recoilless  Rifle       3 

It  is  assumed  that  the  unit  is  engaged  in  a  large  scale  battle 
as  part  of  a  large  force.   The  performance  random  vector  for 
this  unit,  X(t) ,  has  an  outcome  space  ft  ,  consisting  of  the 
eight  binary  3-vectors.   Let  the  state  function  be  defined  as 
follows : 
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s{(l,l,l) }  =  1  s{(l,0,0) }  =  5 

s{ (1,0,1)  }  =  2  s{  (0,1,0)  }  =  6 

s{ (0,1,1) }  =  3  s{ (0,0,1) }  =  7 

s{(l,l,0)  }  =  4  s{  (0,0,0)  }  =  8 

The  decision-maker  is  now  faced  with  the  problem  of 
constructing  the  absolute  effectiveness  function  g(x).   Sup- 
pose the  hypothetical  team  is  considered  to  have  two  missions. 
At  any  time  the  weapon  is  functioning,  the  primary  mission  is 
to  detect  enemy  tanks  and  bring  them  under  fire.   If  the 
weapon  is  in  the  failed  condition,  the  secondary  mission  is  to 
detect  enemy  tanks  and  direct  the  fire  of  adjacent  teams 
against  them.   From  experience,  field  test,  or  some  other 
information  source,  the  user  determines  that  when  both  person- 
nel and  the  weapon  are  intact  the  team  is  capable  of  bringing 
effective  fire  on  eight  tanks,  per  hour.   He  decides  to  base 
his  absolute  effectiveness  function  on  this  fire  capability, 
and  his  first  step  is  to  set  g{ (1,1,1)}  =  8.   He  feels  that 
in  most  units  the  leader  and  gunner  will  have  equal  ability  in 
firing  the  weapon.   Therefore  when  either  member  becomes  a 
casualty  when  the  weapon  is  functioning,  the  decrease  in  effec- 
tiveness is  about  the  same.   The  user  estimates  the  reduced 
capability  to  be  about  five  tanks  per  hour,  and  defines 
g{  (0,1,1)}  =  g{ (1,0,1)}  =  5.   When  the  weapon  fails,  the  team 
performs  the  secondary  fire  direction  mission.   The  user  deter- 
mines that  both  members  together  can  detect  and  bring  fire  on 
as  many  as  three  tanks  per  hour.   However,  if  the  leader  alone 
survives,  he  has  a  capability  of  two  tanks  per  hour,  while  the 
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gunner  alone  has  a  capability  of  only  one  per  hour  because  he 
has  less  experience  and  training  than  the  leader  in  directing 
fire.   In  the  two  cases  where  both  members  are  casualties  the 
capability  is,  of  course,  zero.   Thus  the  decision-maker 
defines  his  absolute  effectiveness  as  follows: 

■  •       g( 

;i,0,l)  }  =  d(2)  =  g5  =  5 


g{ 
g{ 
g( 
g{ 
g( 
g( 
g{ 


1,1,1) }  =  d(i)  =  g6  = 


0,1,1) }  =  d(3)  =  g5  =  5 
1,1,0)  }  =  d(4)  =  g4  =  3 


!1,0,0) }  =  d(5)  =  g3  =  2 

0,1,0) }  =  d(6)  =  g2  =  1 

;0,0,1) }  =  d(7)  =  q1    =  0 

0,0,0) }  =  d(8)  =  q1   =  0 

The  only  requirement  to  be  met  in  practice  is  that  the 

resulting  function  faithfully  reflect  the  true  values  of  the 

decision-maker.   Additionally,  when  several  units  are  being 

compared,  a  separate  function  must  be  defined  for  each  one, 

and  however  they  are  constructed,  it  must  be  the  case  that 

states  of  different  units  adjudged  by  the  user  to  be  equally 

effective  have  the  same  value  of  g  (x)  . 

As  indicated  in  (3-4),  g(x)  can  be  expressed  as  follows: 

g(x)  =  (1-0)  01(x)    +    (2-1)  02(x)  +  (3-2)  p^  (x)  +  (5-3)  jz$4  (x)  + 

(8-5)  0c(x)    where  the  structure  functions  0.(x)  are  defined 
_>  —  l  — 

in  the  table  below. 
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X 

01(x) 

*2<x> 

*3(x) 

04(x) 

*5(x) 

(1,1,1) 

1 

1 

1 

1 

1 

(1,0,1) 

1 

1 

1 

1 

0 

(0,1,1) 

1 

1 

1 

1 

0 

(1,1,0) 

1 

1 

1 

0 

0 

(1,0,0) 

1 

1 

0 

0 

0 

(0,1,0) 

1 

0 

0 

0 

0 

(0,0,1) 

0 

0 

0 

0 

0 

(0,0,0) 

0 

0 

0 

0 

0 

The  block  diagrams  of  coherent  systems  associated  with  these 
structure  functions  are  shown  in  Figure  6. 

The  decision-maker  also  wishes  to  define  a  time  variance 
function  for  the  tank  killer  team  which  enhances  the  actual 
effectiveness  of  every  state  the  longer  the  engagement  lasts. 
He  feels  that  enemy  losses  and  confusion  make  friendly  fire- 
power more  valuable  later  in  the  battle.   He  consequently 
defines 

f(t)  =  2-e_t. 
Thus  f(0)  =  1  and  at  the  outset  h(x,0)  =  g(x),  and  after  one 
hour  h(x,l)  =  (1.633)g(x). 

Assume  the  following  conditional  distribution  for  S(t)  is 

given . 

-7t  /xx     -3t   -4t 

Plx(t)  =  e  Pl5(t)  =  e    "e 

-6t   -7t  ,.,     -2t   -3t 

p12(t)  =  e    -e  Pl6'(t)  =  e    ~e 

,.»     -5t   -6t  ...  -t   -2t 

P13(t)  =  e    -e  p17(t)  =  e   "e 

P14(t)  =  e"4t-e~5t        Pl8(t)  =  1"e"t 
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Figure  6.   Coherent  Systems  Associated  with  the  Absolute 
Effectiveness  Function  of  Example  1. 
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By    (3-8) , 

I(1't)  "    5^ —  e  -5-  e 

to    -m  5  1         -6t  3      -7t  1         -8t 

I(2't)  "   168    "   —  G  +   — *  "   -8~  e 

m  to  -     3  2     o~5t  ,     1       -6t         1       -7t 

T/d    H  .      1  1        -4t  3        -5t  1        -6t 

I(4't)  ~  13 r"  e  +  -3-  e  '  -T"  e 

1(5, t)  =  -2 L.  e-3t  3        -4t   _   _1_     -5t 

l3,t'         60  3  T^  3^  e 

1  -2t  -3t  1        -4t 


1(6, t)    =   -4 e    *u   +   e 


e 


5  _    -t  3         -2t  1         -3t 


e 


1(7, t)    =   -^_  -    2e        +  -J-  e 

1(8, t)    =    2t 1-  +    3e_t ~  e~2t    . 

And  by    (3-7)  , 

8 
q(t)    =     ]T     d(i)     I(i,t) 
i=l 

=81 (l,t)    +51 (2,t)    +    5    1(3, t) 

+  3    I  (4,t)    +    2    1(5, t)    +    1    1(6, t) 

+0    1(7, t)    +    0    I (8,t) 

653  3         -8t        _6_      -7t        _1_      -6t 

280  8e  7e  3e 

_3_      -5t        _1_      -4t        _1_      -3t  -2t 

5      e  4e  3e  e. 

Now  the  given  interval  probabilities  are  equivalent  to 
the  following  reliability  functions  of  the  associated  systems 
in  Figure  6: 
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R^t)  =  e"7t 

R2(t)  =  e~5t 

R3(t)  =  e~4t 

R4(t)  =  e"3t 

R5(t)  =  e"2t 


By  (3-10)  , 


TM  *-\            9      2   Q~7t  x   !    ~8t 
J(l,t)  -  ^ —  e    +  —  e 


T,0  ..7      2    -5t  .   1    -6t 
J(2,t)  =  30-  -  -3-  e     +  -g-  e 


(3,t)  =  j^—  -  — 2~  e     +  — g—  e 


-xiA    *-\  5      2    "3t  .   1    -4t 

J(4,t)  =  jj-  -   -j-  e     +  —q-   e 


J(5,t)  =  -i e  2t  +  -i-  e  3t' 


And  by  (3-9)  , 

5 

q(t)  =   ^   ^j  +  l'^j)  J(J't) 
j  =  l 

=   3  J(l,t)  +  2  J(2,t)  +  1  J(3,t) 
+  1  J(4,t)  +  1  J(5,t) 

653     3    -8t     6    -7t  L   1    -5t 
280     8  7  3 

3    -5t     1    -4t     1    -3t    -2t 
-  — -  e    -  -J-  e    -  -3-  e    -e 

which  agrees  with  the  previously  computed  result. | 
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IV.   UNIT  PERFORMANCE  MODEL 

The  quality  function,  q(t),  has  already  been  introduced 
as  a  measure  of  effectiveness.   The  next  step  is  the  introduc- 
tion of  a  model  to  translate  the  user's  judgment  and  intuition 
about  combat  dynamics  into  a  description  of  the  stochastic 
performance  process  of  the  proposed  unit  structure.   The  pri- 
mary requirements  for  the  model  are  that  it  be  sensitive  to 
structural  alternatives,  and  that  it  provide  results  in  a 
form  suitable  for  input  to  the  quality  function.   There  are 
two  alternative  forms  of  performance  information  which  can  be 
used  in  calculation  of  q(t).   If  expression  (3-9)  is  used,  the 
inputs  required  are  the  reliability  functions  of  the  associ- 
ated coherent  systems.   If  the  calculation  is  performed  accord- 
ing to  expression  (3-7),  the  information  needed  is  the  set  of 
interval  probability  functions  {p. , (t) , . . . ,p,N (t) } .   This 
second  type  of  information  is  the  type  of  performance  descrip- 
tion which  is  most  easily  determined  employing  the  model 
proposed  in  this  section. 

The  unit  to  be  modeled  is  assumed  to  be  a  complex  struc- 
ture which  includes,  as  components,  personnel,  equipment,  and 
channels  of  information.   The  assumption  will  be  made  that 
each  of  the  components  is  either  functioning  or  failed.   At 
the  opening  of  an  engagement  all  components  are  considered  to 
be  functioning.   Some  components  such  as  personnel  are  non- 
renewable and  do  not  return  to  the  functioning  condition  once 
they  leave  it.   Other  components  are  renewable  and  can  return 
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to  the  functioning  condition  repeatedly.   For  instance  the 
information  channel  between  a  squad  leader  and  a  team  leader 
which  consists  of  voice  communication  can  be  interrupted  by 
battle  noise,  then  restored  when  the  noise  subsides.   There- 
fore an  engagement  can  be  looked  upon  as  a  series  of  events 
in  time  which  change  the  state  of  the  unit,  through  attrition, 
interruption,  restoration  or  some  combination  of  these  effects, 
The  state  of  the  unit  remains  unaltered  between  these  events. 
One  approach  commonly  used  to  portray  this  process  is  to 
consider  the  performance  of  each  component  to  be  probabilis- 
tically independent.   This  assumption  is  not  acceptable  for 
this  application  since  it  severely  limits  the  consideration 
of  dependence  in  the  attrition  process. 

In  order  to  integrate  structural  sensitivity  into  the 
model,  it  is  necessary  that  the  interevent  times  and  the  na- 
ture of  each  event  be  considered  to  obey  random  distributions 
which  depend  upon  the  surviving  structure.   The  following 
proposal  for  accomplishing  this  is  derived  in  part  from  the 
work  of  Marshall  and  Olkin  [Ref .  9] . 

The  engaged  unit  is  subject  to  a  hostile  environment  which 
consists  of  the  delivery  over  time  of  one  or  more  types  of 
lethal  projectiles,  including  nuclear  weapons,  artillery,  and 
small  arms  fire.   Two  assumptions  are  made  about  these  fires. 

1)   The  delivery  of  each  type  of  projectile  constitutes 

an  independent  time  homogeneous  Poisson  process, 

having  interarrival  times  distributed  exponentially 

with  expected  value  — -,  where  k  is  the  index  of  the 

k 
weapon  type. 

39 


2)  For  each  unit  state  and  each  weapon  type  the  damage 
caused  by  each  round  obeys  a  probability  mass  func- 
tion:  Pk(j|i)  =  P[Projectile  of  type  k  causes 

transition  to  state  j    unit  is 
in  state  i] . 
From  these  assumptions  it  follows  that  when  the  unit  is  in 
state  i,  the  distribution  of  arrivals  of  projectiles  of  type 
k,  which  send  the  unit  into  state  j,  is  Poisson  having  expo- 
nential interarrival  times  with  expected  value  - — r-r-f- . 

Now  in  general  the  unit  in  state  i  can  be  subject  to 
transition  into  state  j  as  a  result  of  several  threats.   If 
there  are  z  types  of  projectiles  which  can  cause  the  transi- 
tion, then  the  distribution  of  all  such  events  is  the  super- 
position of  z  independent  Poisson  distributions,  which  is 
itself  Poisson,  having  exponential  interarrival  times  with 

1 

mean  — ; ■  i •  i  ■  ■  -i ; ; re ,  ■  I  ■  \ . 

PX(D  !i)^1  +  ...  +  Pz(j  |i)T2 

In  addition  to  attrition  events  there  are  also  interruption 
and  restoration  events  which  describe  alterations  in  condition 
of  the  renewable  components.   It  is  assumed  that  for  each  state 
i  all  functioning  renewable  components  have  an  exponential  time 

to  failure  with  mean  ■= ,  where  j  is  the  state  which  the  com- 

ponent ' s  failure  causes  the  unit  to  enter.   Likewise,  if  state 

i  has  a  failed  renewable  component,  then  ■* —  is  the  mean  of  the 

ik 

exponential  time  to  restoration  and  consequent  transition  into 
state  k. 
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Let  T. .  be  defined  as  the  time  interval  from  the  unit's 
entry  into  state  i  until  an  event  occurs  causing  transition 
into  state  j.   T . .  is  then  the  time  until  an  attrition  event, 
or  an  interruption-restoration  event  of  the  proper  type 
occurs.   This  is  again  a  superposition  of  independent  Poisson 
processes  resulting  in  a  Poisson  distribution  having  exponen- 
tial interarrival  times  with  expected  value, 

E[Tij]  =  P,  (j  li)  t,  +  ...  +  p  (j  I  i)  t   +  6.  .  =  r~ ' 
J  1  J  '    1  z  J  '    z    13     lj 

and  survival  function 


P[T. .  >  t]  =  e    1: 


-(Xj4)t 

>  tl  =  e 
ID 

In  general  the  unit  in  state  i  faces  threats  which  can 
cause  transition  to  several  different  states.   The  outcome  is 
determined  by  the  type  of  the  earliest  state-changing  event  to 
occur.   Let  T   by  defined  as  the  length  of  the  time  interval 
between  transition  a-1  and  transition  a,  so  that 

Ta  =  min{Tll,T12 TiN}, 

if  state  i  was  entered  at  transition  a-1.   If  transition  from 
state  i  to  state  i  is  not  possible,  then  T. .  is  considered  to 
be  infinite,  with  X. .  =  0.   Howard  [Ref.  6]  refers  to  this 
type  of  model  as  a  "competitive  process"  since  it  is  equivalent 
to  a  process  in  which,  upon  entry  into  state  i,  the  random 
variables  T.  .  are  sampled  and',  in  effect,  compete  for  the 
earliest  outcome. 
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The  survival  function  of  T   is  derived  as  follows: 

Ut 

P[Ta    >    t]     =    P[Ti;L    >    t,     T.2    >    t T.N    >    t] 

=    P[Ti;L    >    t]        P[Ti2    >    t]...P[TiN    >    t] 

"(Xil    +   Xi2    +    '••    +    AiN}t 
=   e 

Thus  Ta  is  distributed  exponentially  with  mean, 

Thus  the  original  assumptions  of  Poisson  arrival  of  projectiles, 
and  exponential  interruption-restoration  times,  imply  that  the 
times  between  state-changing  events  are  also  exponentially 

distributed  with  mean  r —  . 

i 

The  validity  of  these  assumptions  must  be  considered  in 
any  application  of  the  model.   The  assumption  of  Poisson  inci- 
dence of  fires  is  characteristic  of  relative  stability  in 
delivery  rates.   This  distribution  would  not  be  consistent 
with  scenarios  which  included  large-scale  enemy  reinforcement 
or  temporary  massing  of  fires. 

The  exponential  distribution  of  the  interevent  times,  T  , 
implies  the  following  result: 

-UjHt+s)    _x.s 

P[Ta  >   t  +  s|T0  >  t)  =  2-^ =  e   1 

e 

=  P[T   >  s]  . 

a 

This  result  is  referred  to  as  the  "memoryless"  property  or,  in 
reliability  theory,  as  the  "no  wear"  property.   The  effect  of 
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the  property  in  this  model  is  that  a  unit's  probability  of 
successfully  remaining  in  its  current  state  for  one  addi- 
tional minute,  is  completely  independent  of  the  amount  of 
time  already  spent  in  the  state.   This  property  fails  to  ac- 
count for  such  factors  as  fatigue  or  depletion  of  supplies. 
However  in  any  case  where  the  effect  of  these  factors  is 
insignificant  compared  with  the  actual  attrition  and  interrup- 
tion, the  "no  wear"  property  may  prove  to  be  acceptable. 

The  major  consequence  of  the  exponential  assumptions  is 
that  the  Markov  property  holds.   That  is,  for  t  >  u  >,  0 , 

P[S(t)  =i|s(v),  u>.v>lO]  =P[S(t)  =  i|S(u)  ]  , 

or  the  distribution  of  future  S(t)  depends  only  upon  the  cur- 
rent state,  and  not  at  all  upon  earlier  information.   The 
disadvantage  of  this  property  is  that  it  does  not  permit  con- 
sideration of  such  phenomena  as  momentum  and  collapse  of 
situation.   In  applications  where  such  effects  are  considered 
important,  alteration  would  be  required. 

The  fortunate  consequence  of  the  Markov  property  is  that 
it  qualifies  the  process,  {S(t),  t  :>  0},  as  a  continuous  para- 
meter Markov  chain,  a  type  of  process  which  possesses  a  well 
developed  method  of  solution. 

The  decision-maker  using  this  model  is  presumed  to  have  in 
mind  a  specific  unit  structure  and  a  specific  engagement 
scenario.   From  the  scenario  he  determines  the  variety  of 
projectile  types  to  which  the  unit  is  vulnerable  along  with 
rates  of  fire  t,  .   Then  for  each  state  i  he  must  specify  the 
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mass  functions  P,  (j|i)  for  each  type  threat  k.   These  mass 
functions  are  the  primary  ingredient  of  the  model's  sensi- 
tivity to  structure.   They  should  be  based  on  such  variables 
as  each  individual's  level  of  training,  the  capacity  in  which 
he  is  serving  due  to  attrition  of  key  personnel,  the  degree 
to  which  he  must  expose  himself  to  hostile  fire  as  a  result 
of  his  current  position,  and  the  current  condition  of  channels 
of  communication.   The  specification  of  the  mass  functions 
will  be  the  most  difficult  part  of  parameterizing  the  model. 
But  it  is  indispensible  if  worthwhile  investigation  of  unit 
structure  is  desired. 

The  fact  that  the  Markov  property  holds  leads  to  the 
derivation  of  two  sets  of  simultaneous  first  order  differen- 
tial equations  which  can  be  solved  for  all  the  functions 
p.  .  (t)  .   These  two  sets  of  equations  are  known  as  Kolmogorov's 
forward  equations  and  backward  equations.   The  actual  deriva- 
tion is  not  presented  here  but  is  given  in  Ross  [Ref.  11]. 
Before  presenting  these  equations,  several  notational  concepts 

are  introduced. 

2 
There  are  N   distinct  functions  p. . (t) ,  which  appear  in 

the  functional  matrix 


P(t)  = 


pn(t) 


Pni 


(t) 


PlN(t) 


PNN(t) 
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where  the  function  p. . (t)  appears  in  the  ith  row  and  jth 
column.   For  any  t  >.  0  the  rows  of  p(t)  add  to  one.   Since 

p±j(0)  =  1    i-j 
=0     i^j, 

it  follows  that 


P(0)  =  I  = 


1  0 
0  1 


Lo  . 


.  o 


1  o 
o  lj 


The  functional  matrix  P'(t)  is  defined  as  follows 


P' (t)  = 


P'll(t)  •••  P'lN(t) 


P'Nl(t)  ■•■  P'NN(t) 


where 


p' • • (t)  = 
r  l] 


d  p. . (t) 
dt 


Using  the  parameters  already  defined  in  this  section,  the 
constant  matrix  A  is  defined  as  follows: 


A  = 


X 


21 


12 

0 


.  X 


IN 


•Nl 


whe 


re  X .   is  the  reciprocal  of  the  mean  of  T. ..  and  X . .  appears 
ith  row  and  jth  column.   Since  no  transitions  are  considered 
from  any  state  into  itself  X..  =  0,  for  i  =  1,...,  N. 
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The  constant  matrix  B  is  defined  as  follows: 


X12  * ' '  X1N 


B   = 


X21  "  X2 


Nl 


-X 


N 


N 


where  X .    =  )  X  .  . ,    the  reciprocal  mean  of  the  interevent  time 
1    Z_.  13  r 

j=l     when  the  unit  is  in  state  i. 
Kolmogorov's  equations  are  stated  in  matrix  form  below. 
The  backward  equations  are  given  by 


(4-1) 


P*  (t)  =  B  P(t)  , 


with  initial  condition  P(0)  =  I.   The  forward  equations  are 
given  by 


(4-2) 


P'  (t)  =  P(t)  B, 


with  initial  condition  P(0)  =  I.   The  operation  indicated  on 
the  right  of  each  expression  is  matrix  multiplication.   The 
fact  that  B  P(t)  =  P(t)B  is  of  special  interest  since  matrix 
multiplication  is  not,  in  general,  commutative.   The  solution 
of  the  process  will  consist  of  solving  the  backward  equations 
using  Laplace  transforms. 

The  Laplace  transform  of  any  functional  matrix  A(t),  is 
indicated  by 

L{A(t) }  =  A(s)  , 

which  is  simply  the  matrix  of  Laplace  transforms  displayed  in 
the  same  order  as  in  A(t).   For  any  functional  matrix  A(t), 
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(4-3)         L{A'(t)}  =  A'(s)  =  sA(s)  -  A(0), 

an  equation  proved  by  Kreyszig  [Ref .  8] .   Also  for  any  func- 
tional matrix  A(t)  and  any  constant  matrix  C, 

(4-4)  L{C  A(t) }  =  C  A(s) 

(4-5)  L{A(t)C}  =  A(s)C 

If  Laplace  transforms  are  taken  of  both  sides  of  (4-1) 
the  result  is 

L{P'  (t)  }  =  L{B  P(t)  }  , 
which  by  (4-3)  and  (4-4)  is  equivalent  to 

sP(t)  -  P(0)  =  B.  P(t)  , 

or  since  P(0)  =  I, 

sP(t)  -  B  P(t)  -  I, 
or 

[si  -  B]  P(t)  =  I. 

Premultiplying  both  sides  by  the  inverse  of  [si  -  B]  gives, 

P(t)  =  [si  -  B]"1, 

which  gives  the  solution  for  L{P(t)}  in  terms  of  the  constant 
matrix  B.   The  solution  is  complete  when  the  transform  is 
inverted  to  give  P(t). 

Howard  [Ref.  6]  shows  that 


if  and  only  if 


L{P(t) }  =  [si  -  B]"1, 


P(t)  =  eBt  , 
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where 

2  3 

_Bt         T         ^  B       ,2     .    B         3 

e=I+Bt+yTt      +   ~—  t      +...        , 

is  called  the  matrix  exponential  function.   This  expression 
may  be  used  to  calculate  approximations  to  the  elements  of 
P(t),  but  will  not  provide  exact  solutions.   Such  a  solution 
can  be  arrived  at  by  finding  the  inverse  of  [si  -  B]  and 
using  the  method  of  partial  fractions  to  translate  the  result 
into  invertible  Laplace  transforms. 
Example  2 

Consider  again  the  tank  killer  team  of  Example  1,  with  the 
same  three  components,  all  considered  nonrenewable,  and  the 
same  state  function  s (x) . 

The  scenario  of  interest  involves  two  types  of  external 
threat.   The  first  threat  consists  of  small  arms  projectiles, 
delivered  with  a  Poisson  rate  of  fire  t,  =  50/hour.   The 
second  threat  is  artillery  fire  with  rate  t2  =  10/hour. 

The  decision-maker  now  determines  the  approximate  proba- 
bility mass  function  for  each  state-threat  combination. 

For  state  one  and  the  small  arms  threat,  the  user  first 
notes  a  single  round  cannot  cause  more  than  one  casualty.   He 
therefore  defines  P1(5|l)  =  P1  (6  |  1)  =  Px  (7  1 1)  =  P]L  (8  1 1)  =  0, 
since  these  transitions  would  entail  multiple  attrition.   He 
feels  that  the  gunner  is  about  50%  more  vulnerable  to  small 
arms  than  the  leader  because  of  his  firing  position,  and  the 
fact  that  the  backblast  of  his  weapon  pinpoints  his  location. 
The  weapon  is  considered  only  half  as  vulnerable  as  the  leader 
because  its  small  presented  profile.   The  user  thus  estimates 
the  following  mass  function: 
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P1(l|l)  =  .94  P1(5|l)  =  0 

Px(2  |1)  =  .03  P1(6  |1)  =  0 

P1(3|l)  =  .02  P1(7|l)  =  0 

P1(4|l)  =  .01  P1(8|l)  =  0 

where  the  actual  magnitudes  may  be  based  on  field  data, 
simulation,  or  intuition. 

For  state  one  and  the  artillery  threat  multiple  attrition 
is  possible.   The  user  feels  that  each  artillery  round  has  a 
capability  to  destroy  all  three  components  and  defines 
P„  (8  1 1 )  =  .1  .    The  gunner  is  considered  more  vulnerable  to 
artillery  fire  for  identical  reasons  as  cited  under  small  arms 
The  decision-maker  defines  this  mass  function  as  follows: 

P2(l|l)  =  .60  P2(5|l)  =  .1 

P2(2|l)  =  .05  P2 (6  1 1)  =  0 

P2(3|l)  =  .1  P2(7|l)  =  0 

P2(4  |1)  =  .05  P2(8  |1)  =  .1   . 

The  probability  of  transition  into  state  seven  is  considered 
to  be  zero,  since  a  round  placed  well  enough  to  eliminate  both 
personnel  would  be  expected  to  destroy  the  weapon  as  well. 
Transition  to  state  six  is  similarly  considered  impossible. 
The  actual  magnitudes  are  again  based  on  the  user's  judgment. 

Since  there  are  no  renewable  components  the  parameters  A  . . 
can  now  be  computed : 
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A,,  =  0     (trivial  transitions  are  ignored) 

X12    =  P1(2I1)T1    +    P2(2|Dt2 

=  (.03)  (50)    +    ( .05)  (10) 

=  2 


X13    =    Pi^3!1)1!    +    P2(3|Dt2 


=    (.02)  (50)    +    (.1)  (10) 
=    2 


X14    =   P1(4|1)t1    +   P2(4|1)t2 


=    (.01)  (50)    +    (.05)  (10) 
=    1 


X15    =    P1(5|Dt1    +   P2(5|1)t2 


=  (0)  (50)    +    (.1)  (10) 

=  1 

A16    =  Pl(6 '1)T1    +    p2(6I1)t2 

=  0 


X17    =    Pl(7 l1)Tl    +   P2(7'1)T2 


=    0 


X18    =    P1(8|Dt1    +    P2(8|1)t2 


=    0    +     (.1) (10) 
=    1 
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The  construction  of  the  full  set  of  mass  functions 
results  in  the  constant  matrix  A  shown  in  Figure  7,  and  the 
functional  matrix  [si  -  B]  in  Figure  8. 

The  method  of  solution  is  to  invert  [si  -  B]  by  the  method 
of  cofactors  using  the  relation 

(-l)i+j 

p..(s)  =  4—^ 5-rM..(s)  . 

*i;j       [si  -  B  |   31 

where  |sl  -  b|  is  the  determinant  of  the  matrix  [si  -  B]  ,  and 
M.. (s)  is  the  determinant  derived  from  Isl  -  b|  by  deleting 

jl  \  \         j  r> 

the  jth  row  and  ith  column. 

The  advantage  of  the  cofactor  method  is  that  it  permits 
the  elements  of  [si  -  B]    to  be  computed  individually  saving 
time  in  cases  where  only  a  portion  of  the  inverse  is  needed. 

The  results  are  as  follows: 


pll(s)  =  FTT  =  L{e~7t}' 


~  (    ^         2        1      1 
p12ls;    (s+7) (s+5)  "  s  +  5    s+7 

T  f    -5t     -7t, 
=  Lie     -  e    i , 


P13(S}    (s  +  7)  (s+6)  ::  s  +  6    s+T 
=  L{2e~6t  -  2e"7t}' 


~  .    .         1        1      1 
P14(S;    (s+7) (s+6)  "  s+6  "  s+7 

_  r      -6t        -7t, 

=  Lie     -  e    } , 
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Figure  7.   Constant  Matrix  A 
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°j 


s+7  -2  -2  -1  -1 
0  s  +  5  0  0  -2 
0       0      s+6      0      -1 


L0 


s+6  0  -1 
0  s+4  0 
0       0      s+4 


0 

-l" 

0 

-1 

-3 

-2 

-3 

-2 

0 

-4 

0 

-4 

s  +  4 

-4 

0 

s 
1 

Figure  8.   Functional  Matrix  [si  -  B] 
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~  1  4  2 

P15(s)    (s+7) (s  +  4)  +  (s+7) (s+5) (s+4)  +  (s+7) (s+6) (s+4) 


2     2 


s+7    s+6    s+5    s+4 


=  Lie     -  e    -  2e     +  2e    } , 


P16(S)    (s+7) (s+5) (s+4)  +  (s+7) (s+6) (s+4) 

1  3 

1      2      2      2 


s+7    s+6    s+5    s+4 


_  T/  -7t    1   -6t    0  -5t    3   -4t, 
=  L{e     ~   Te  -2e     +  t  e    '  » 


P17(S)    (s  +  7)  (s  +  6)  (s  +  4) 

9_      3 

3      2   +   2 


s+7    s+6    s+4 

T  IO    -It         9   -6t  3   -4t-, 

=Li3e     ~    2   e  2   e     ' 

~   r  .  48 32    6 

p18l  }  (s+7)  (s+6)  (s+4)s  (s+7)  (s+5)  (s+TTs    (s+7)  (s  +  6)s 

+       2       .       4  1 

(s+7) (s+5)s    (s+7) (s+4)s    (s+7)s 

1      2      3  3      5 

+  r-s-     + 


s     s+7    s+6    s+5    s+4 

T/,  ~  -7t  ,  -  -6t  ,  .,  -5t    c  -4ti 
=   Lil--2e     +  3e     +  3e     -  5e    }. 

8 

As  required  y   p,  .(t)  =1,  t  >.  0 . 
j=l 

The  interval  probabilities  {p, , (t) , . . . ,  p.g(t)}  are  now 
known  and  the  quality  function  can  be  computed  for  any 
functions  f(t),  and  g(x)  or  d(i).| 
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It  should  be  noted  that  the  example  was  particularly 
simplified  by  the  fact  that  the  matrix  [si  -  B]  had  no  non- 
zero entries  below  the  diagonal.   The  reason  for  this  is 
that  all  three  components  were  considered  to  be  nonrenewable. 
Thus  when  the  unit  left  any  state  it  could  never  return  to 
it.   And  by  arranging  the  states  in  order  so  that  the  earlier 
states  corresponded  to  the  higher  rows,  the  triangular  form 
was  assured.   The  resulting  denominator  of  all  the  expres- 
sions p. . (s)  was  of  the  form 

8 

|si  -  b|  =  7T(s+x.) 

i=l    1 

which  is  already  in  the  factored  form  necessary  for  the  partial 
fraction  expansion. 

If  however  at  least  one  of  the  components  had  been  renew- 
able, there  would-  have  been  at  least  one  positive  entry  below 
the  diagonal  and  the  denominator  of  the  expressions  p. . (s) 
would  have  been  of  the  form, 

8 
|sl  -  B|  =  ~JT  (s+X  .)   -  K(s) 
i=l 

where  K(s)  is  a  polynomial  function  of  s  of  degree  less  than 
eight.   Therefore  | si  —  B |  is  a  polynomial  of  degree  eight 
whose  roots  and  factors  are  unknown  in  general.   This  computa- 
tional problem  is  formidable  since  the  degree  of  polynomials 
to  be  solved  will  equal  N  =  2   (where  n  is  number  of  system 
components) ,  and  there  is  no  exact  solution  for  polynomial 
equations  of  degree  greater  than  four.   There  are  however 
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computer  routines  which  perform  approximation  procedures, 
which  may  prove  to  be  of  value  in  employing  this  model. 

Finally,  a  rough  measure  of  effectiveness  may  be  computed 
from  q(t)  without  inverting  the  Laplace  transforms  when 
f (t)  =1,  t  >  0 .   Let  the  quality  index  be  defined  as  follows 

JT  00 

f       r 

Q  =  q(»)  =  /_d^)  J    P±j  (u)du 
j-1     o 

00 

d( j)  /  e~   pi .  (u)du 
j  =  l     o 

d(i)   d. 


Thus  Q  is  the  expected  aggregate  value  over  an  infinite  time 
period.   Of  course,  any  realistically  defined  unit  will 
reach  a  state  with  zero  effectiveness  in  finite  time,  so 
that  Q  has  a  finite  value.   However  this  measure  does  not  dif- 
ferentiate between  units  whose  quality  functions  are  thinly 
distributed  over  the  real  time  line  and  those  whose  effec- 
tiveness is  heavily  concentrated  in  short  intervals  of  reason- 
able engagement  duration. 
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V.   CONCLUSION 

The  purpose  of  this  thesis  has  been  to  outline  a  proposed 
modeling  viewpoint.   It  has  not  been  advanced  as  a  finished 
product.   It  is  hoped  that  the  concepts  introduced  will  be 
considered  as  a  worthwhile  point  of  departure  for  additional 
research,  which  could  result  in  establishment  of  a  valuable 
decision  tool  for  choosing  among  alternative  unit  structures. 

The  most  critical  unsolved  problem  is  the  computational 
difficulty  involved  in  solving  the  model.   In  addition  to  the 
inversion  of  Laplace  transforms,  the  current  solution  neces- 
sitates the  evaluation  of  large  order  determinants,  which  can 
become  unwieldy  even  with  electronic  computers.   Thus  the 
current  span  of  application  is  limited  to  small  organizations. 
Research  to  extend  this  span  should  concentrate  on  finding 
simplified  computational  procedures  for  exact  solution,  or 
alternatively,  determining  efficient  and  accurate  approximation 
methods.   The  use  of  the  infinite  series  definition  of  the 
matrix  exponential  function  as  an  approximation  method  was 
introduced  in  Section  IV.   Unfortunately,  this  expression 
requires  sequential  multiplication  of  large  order  matrices,  an 
operation  which  is  as  computationally  explosive  as  evaluation  of 
determinants . 

Current  research  in  reliability  theory  may  offer  a  valuable 
approximation  method,  particularly  since  all  of  the  random 
processes  and  measures  of  effectiveness  were  shown  to  be 
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expressible  as  linear  functions  of  variables  related  to 
reliability  of  coherent  systems.   Esary  and  Proschan  [Ref.  4] 
show  the  usefulness  of  reliability  bounds  and  approximations 
based  on  the  concepts  of  minimal  cut  sets  and  minimal  path 
sets  of  components.   Adaptation  of  this  method  to  the  avail- 
able type  of  data  inputs  may  prove  useful  to  the  current  model. 

A  possible  exact  solution  method  using  flow  graph  analysis 
is  proposed  by  Howard  [Ref.  6] ,  but  at  present  it  appears  to  be 
subject  to  an  even  greater  computational  burden  than  the  cur- 
rent method. 

In  order  to  enhance  the  accuracy  of  the  model  it  may  be 
desirable  to  alter  the  assumption  of  exponentially  distributed 
interevent  times  in  favor  of  some  other  distribution.  This 
would  change  the  process  to  one  of  the  semi -Markov  variety. 
These  processes  are  more  difficult  to  solve  in  general,  but 
offer  greater  flexibility  in  definition.  Ross  [Ref.  11]  and 
Howard  [Ref.  6]  consider  several  solution  methods. 

A  highly  desirable  addition  to  the  model  would  be  a  provi- 
sion for  random  threat.  Rather  than  assuming  constant  Poisson 
rates  of  enemy  fire,  such  a  refinement  would  allow  rates  to 
vary  randomly,  thus  representing  periods  of  non-engagement  as 
well  as  change  of  situation  during  engagements.  The  analogous 
reliability  concepts  of  random  wear  and  random  shock  have  been 
studied  by  Gaver  [Ref.  5]  and  Reynolds  and  Savage  [Ref.  10]. 

Despite  its  capacity  for  refinement,  the  proposed  model  in 
present  form  does  fill  many  of  the  needs  for  the  study  of  struc- 
ture.  It  offers  the  decision-maker  a  model  which  uses  his  own 
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judgment,  in  an  intuitively  clear  process  of  parameterization, 
in  order  to  compute  a  functional  measure  of  effectiveness 
based  on  fundamental  functions  which  he,  himself,  defines. 
Indeed,  with  the  exception  of  clearly  stated  assumptions  made 
in  building  the  model,  and  the  purely  mathematical  solution 
process,  the  model  in  use  is  completely  established  by  the 
user,  with  whom  ultimate  decision  responsibility  rests. 
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